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Abstract — Renewable energy sources can be attached to sensor 
nodes to substantially improve the performance of sensor net- 
works. In networks with renewable energy sources, conservative 
energy expenditure may lead to missed recharging opportunities 
due to the batteries being full, while aggressive usage of energy 
may cause the network to be intermittently disconnected and 
reduction in coverage. Thus, new techniques and algorithms 
must be developed for sensor networks with replenishment to 
balance these seemingly contradictory goals. In this paper, we 
consider a sensor network with renewable energy sources under 
the SINR interference model, where the achievable rate of each 
link depends on its SINR. We develop a provably efficient and 
distributed solution to maximize the total sensing rate subject to 
quality of service constraints. Our solution is a joint rate control, 
power allocation, and routing algorithm which is computationally 
simple and applicable to non-rechargeable sensor networks as 
well. 

I. Introduction 

Wireless sensor networks have many applications in mon- 
itoring, maintenance, and environment sensing [1], [2], [3]. 
New developments in the area of renewable energy [4], [5] 
will help facilitate earlier and more wide-spread deployment 
of these networks by increasing their operational capacity and 
lifetime. These renewable sources of energy could be attached 
to the nodes, and would typically replenish the batteries at 
a slow rate (e.g., compared to the rate at which energy is 
consumed by a continuous stream of packet transmissions) 
that could be variable and dependent on the surroundings. 

Unlike traditional battery-limited sensor networks, in net- 
works with replenishment, the network lifetime can be infinity, 
and a node with a full battery may not be necessarily desirable. 
This is because a node with a full battery implies missed 
recharging opportunities, which in turn implies a reduced 
feasible sensing and transmission rates. On the other hand, an 
over aggressive use of energy may lead to lack of coverage or 
connectivity for certain time periods, or result in the network 
being incapable of transferring time-sensitive data to the sink 
because certain nodes have depleted their batteries. Thus, new 
techniques and algorithms must be developed for networks 
with replenishment to balance these seemingly contradictory 
goals. In this paper, we develop a distributed joint rate control, 
power allocation and routing algorithm, which is provably 
efficient for rechargeable networks. We use the so-called 
"physical interference model," for which the transmission rate 
is a function of the signal to interference and noise ratio 



(SINR) and influenced by the power vector over the entire 
network. 

While the problem of energy management in wireless net- 
works has seen considerable attention, there have been few 
works [6], [7], [8], [9], [10], [11] that also include energy 
replenishment. In [6], [7], [8], [9] the authors study the 
problem of dynamic node activation, energy harvesting, sleep- 
wake scheduling, and duty cycle controlling respectively under 
rechargeable networks. In [10] and [11], motivated by [12], 
the authors consider jointly managing the data and battery 
buffers to maximize the throughput with renewable energy 
sources under protocol interference model. The coupling be- 
tween the data and battery buffers is what makes the problem 
notoriously difficult to solve using standard optimization based 
approaches. For example, unlike [12], which utilizes the fact 
that a static allocation policy is optimal, it is not even clear 
whether an optimal policy would be static in a rechargeable 
network setting. Specifically, with a battery buffer, there is an 
additional energy constraint that the allocated energy should 
be within the battery state, and this constraint is even more 
difficult than the average power constraint. In [10], the authors 
consider infinite data buffer and finite battery buffer sizes. 
They assume that the replenishing process is i.i.d, and show 
that the probability of battery state being less than the peak 
power or close to the full battery state vanishes as the battery 
size grows with their algorithm. 

All the prior work on networks with replenishing nodes rely 
on discrete interference models, in order to develop practical 
solutions to the optimal resource allocation problem. Despite 
their simplicity, discrete interference models could be poor 
approximations to the actual interference, observed in wireless 
networks. On the other hand, physical interference models 
such as the SINR model are known to more accurately capture 
the impact of the interference more closely. Prior work on 
resource allocation and utility optimization under the SINR 
model include [13], [14], [15], [16], [17], [18]. In [13], [14], 
[15], the authors formulate optimization problems with QoS 
constraints on SINR. In [16], [17], [18], the authors consider 
an optimal power allocation problem such that a utility func- 
tion of SINR is maximized. They have assumptions such as 
high or low interference regime to avoid the non-convexity. 
None of these studies consider energy replenishment. In this 
paper, we consider the power allocation problem combined 
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with rate control, to achieve the maximal sensing rate with 
energy replenishment using the SINR model. We do not 
assume high or low interference when use the rate power 
function, and our solution is applicable to non-rechargeable 
wireless networks as well. 

The main contributions of the paper can be summarized as 
follows: 

• We develop a mathematical framework to handle the 
coupling between the data and battery queues that exists 
in networks with replenishment. Our mathematical frame- 
work allows for multihop, dynamic routing, and does not 
require future knowledge of the data and replenishment. 

• Unlike existing works on utility maximization and net- 
work optimization, in our framework, we allow for a 
general arrival and replenishment process that does not 
require stationarity or ergodicity. This is especially impor- 
tant in such networks, because the replenishment profiles 
are typically non-stationary. 

• We develop a fully distributed, computationally simple, 
algorithm that achieves a provable constant fraction of 
the optimum solution under the SINR interference model. 
Note that the development of provably efficient dis- 
tributed algorithms under the SINR interference model 
has been an open problem even for networks without 
replenishment. 

II. System Model and Problem Formulation 
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Fig. 1. Multihop Network Model 

We consider a rechargeable wireless sensor network with 
TV nodes and L links, as illustrated in Figure. 1. Each node 
n G Ji = {1,2, ...,7V} is attached to a renewable power 
source for replenishment. We assume a time slotted system. 
In time slot t, A n , e (t) and R n ^ e (t) denote the amount of 
available data for sensing and the actual amount of sensed 
data, at node n that are destined to node e. Hence, R n , e (t) < 
A n ,e{t), Vt > 0, Vn,e e K. Data flows may traverse 
multiple hops to the destination and their routes are not fixed. 
Let r n (t) denote the amount of replenishment at node n 
in slot t. We assume that limsupy^^ ^ 2~2t=o A n ,e(t) < 

A max , HmSUPj^ yEL'o 1 A lA f ) ^ A max, (0 < A max < 

oo, 0 < A 2 m . ix < oo) and 0 < r n (t) < r max , (0 < 

^max ^ 

oo), Vt > 0, Vn, e€N, Without loss of generality, we assume 
the instantaneous amount of data A n _ e (t) and replenishment 
r n (t) are available at the end of slot t. We assume that 
each node n maintains a separate infinite data buffer with 



state q% e (i) for the flow destined to e, and also maintains 
a finite battery buffer with size B h n < oo and state q^(t). Note 
that we focus on infinite data buffer size and finite battery 
buffer size for ease of illustration. Following our approach in 
[11], we can also allow for finite battery and data buffers. 
We assume that, for each link I e L = {1,2, ...,L}, the 
power used to transmit a packet by the transmitting node 
and the power used to receive and decode a packet by a 
receiving node are identical and equal to Pi{t). Note that 
this assumption is merely for notational convenience and 
the generalization to the asymmetric case is straightforward. 
Vector P(t) = {P 1 (t),P 2 {t),...,P L (t)} denotes the power 
assignments over all the links at time t. With this power 
allocation, the achieved data rate at link I is m(P(t)) in 
that time slot. Let P be the sequence of power assignments 
{P(0), P(l), . . . , P(T - 1), . . .} and R be the sequence of 
sensing rates {R(0),R(1), . . .,R(T-1), . . .}. Let fi„ denote 
the set of directed links that originate from node n, and let 
8„ denote the set of directed links that terminate at node n. 
We assume that the total power used by a node in each slot is 
upper bounded by the peak power -P pea k of sensor hardware, 
i- e -' 2~2ien n u e„ ^K*) — -Ppeak- The rate-power function is 
assumed to have the following standard form under the SINR 
interference model: 



W (P(t))=Iog 1 + 



fftran(0,rec(/)-P«(*) 



log 1 



^rec(0 + YsieCijtl 9txsa(i),rec(l)Pi{t) 

W) \ 



(1) 



where gtran(i).rec(j) is the channel gain from the transmitting 
node tran(i) of link i, to the receiving node rec(j) of link j, 
N n is the background noise at node n, o\ = '— — is the 

° 9tran(0,rec(!) 

noise to gain ratio of link I. Note that, even though log(-) is a 
concave function, ni(P(t)) is not a concave function of P(t). 
We can see from this rate-power function that the transmission 
rate on link I is influenced by the power assignment on all 
other links over the entire network (interference) based on the 
SINR model. 

Definition 1: p° t is the long-term battery discharge ratio for 
node n: 



p° n = limsup 



(2) 



where 

Definition 2: 



ln(t) = indicator that battery is discharged 
in slot t for node n 

0 if£j e n nU e B -fi(*)<9n(*) 



1 otherwise 



(3) 



Definition 1 is a suitable QoS metric for rechargeable sensor 
networks that represents how often the battery is discharged. 
Some applications may require the battery state to be always 
above certain positive level, then Definition 1 can be easily 
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modified to represent how often the battery is below the 
desired level, and our solution structure works as well. Note 
that the ratio p° represents the actual probability, if the system 
variables were to satisfy the strong law of large numbers 
(SLLN), which we do not assume here. From a quality of 
service perspective, it is desirable to keep p° to be less than 
some small upper bound r]°. 

Our overall objective is to maximize the long-term average 
total sensing rate at all nodes destined to all destinations 
subject to the constraints on the "rate" of battery discharge for 
all nodes. With the above setting, we formulate the problem 
as follows: 



T-l 



(A) maxliminf ^ V V RnA 1 ) 

1 ,ft t=0 n.ee'N 

< e (*+ 1) < (< e (*) - E tf(P(t))) + + Rn,e(t) 



+ £>f (£(*)), n^e, 



q b n {t + l) = min 



<&(*)- W) + r n (t),B b n 

zeo„ue„ 



0< <min[^(t),P peak ], 
iefi„ue„ 

N 

J2tf(P(t)) = m(P(t)), R n , e (t) < A n , e (t), 

e=l 

1 T_1 

limsup- ^2qi !e (t) < oo, n ^ e, 



T^oo 



P° n <Vn, VneK. 



(4) 

(5) 
(6) 

(7) 

(8) 
(9) 



In Problem (A), constraints (4) and (5) describe how the data 1 
and battery queues evolve, respectively. Constraints (6) are 
the energy conservation equations ensuring that we cannot 
oversubscribe the energy that is unavailable in the battery, 
nor can we exceed the peak power level. Constraints (7) 
are the rate conservation equations that bound the actual 
amount of sensed data R n ^ e (t) by the available amount of data 
A n , e (t), and share the transmission rate of a link among all 
the destinations in slot t. Constraint (8) is the QoS constraint 
for the data queue: we need to keep all the data queues stable. 
Constraint (9) is the QoS constraint that guarantees a battery 
discharge rate of r]°. 

We proceed along the following lines: 

• We define virtual queues for the battery buffers. We 
show that keeping these virtual queues stable ensures that 
Constraint (9) is met. 

• We design a distributed algorithm that keeps the data 
queues and virtual battery queues stable, and at the same 

'Note that if the physical layer does not allow full duplex operation, g„ t „ 
can be chosen to be infinite, so that the node will not be able to transmit and 
receive at the same time. 



time achieves a constant fraction of the global optimum 
of Problem (A) within a provably small gap. 
Definition 3: We define virtual battery queues Vn G N: 

q b n (t + 1) = ($(t) - V° n ) + + l°(t) r n (t) + M n (t) 

+ E + , (10) 

zefi„ue„ 

where M n (t) = (q b n (t) - E le n n ue n W + r„(t) - B h n ) + 
is the amount of missed replenishment and 7°(t) are defined 
in Equation (3). 



r.<f) 



Fig. 2. Battery Queue and Virtual Battery Queue 

From Figure 2, J2ien n ue n Pl(t)-r n (t)+M n (t) is actually 
the decrement of <j£(t) from slot t to t + 1. Since the battery 
has finite size, this term vanishes on average as time goes 
to infinity. Then, p° = limsupy^^ ^ l - a T n and 77° can 
be viewed as the long-term input and output rate of 
respectively. Thus, it is reasonable to expect that q\(t) being 
stable implies p° < r]°. 

Theorem 1: If all the virtual battery queues <7„(t) are 
strongly stable, i.e., limsupj,^^ ^ J2t=o Qn(*) < 00 > ^ n *= 
INT, we then have p° < 77°, Vn G N. 

The proof can be found in Appendix A. We give the main 
algorithm in the following section. 

III. Joint Rate Control, Power Allocation and 
Routing Algorithm for Multihop Networks 

We begin by defining the notion of crosstalk. 
Definition 4: 



hi. 



5tran(i),rec(j) 
Strang) ,rec(j) 



is the crosstalk between a given link i to another link j. 

From Equation (1), the rate on any link / is a function of the 
power allocation over all the links in the network. To localize 
the power control, our power allocation scheme activates, i.e., 
allocates a non-zero power only to links with crosstalks under 
a certain threshold H > 0. 

Definition 5: The conflict set £; of link I is defined as 
follows: 

£; = {k G £ : k^l; h k ,i > H or h Uk > H} 

We will discuss how we choose H later. If link I is activated, 
the links in its conflict set £j can not be activated. Note that 
for any two links i and j, if i G £j, then j G 

Definition 6: The conflict degree K(l) of link I is the 
maximum number of links in £; that can be activated simul- 
taneously. The conflict degree K of the network is defined as 
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follows: 



K 



maxK(l). 



The notion of conflict degree is similar to the interference 
degree under protocol interference model [19]. We use the 
term "conflict" here because even links that are not in each 
other's conflict list can still interfere with each other under the 
SINR model. 

Definition 7: A set of links S is a maximal schedule if all 
the links are not conflicting with each other, but adding any 
link outside 5* will result in a conflict. 

Note that J2kei\J £ ; ] {kes] > 1, VZ e £. 

Our algorithm is composed of three components: rate con- 
trol, power allocation, and routing. 

Multihop Rate Control (MRC): 
If e (t) < , node n chooses to sense all the avail- 
able data until reach a given bound W, i.e., i? n e (t) = 
min[A„ ;e (t), W\, otherwise, reject all the data, i.e., R n ^(t) = 
0. Here,V, W e (0,oo). 

Multihop Power Allocation (MPA): 
Here the goal is to ensure that no node transfers data of a 
flow to a relay node that is not the destination of that flow, 
unless the differential backlog for that flow is greater than 
a fixed value 7 > 0. We will choose the value of 7 such 
that the resulting backlog of the receiving node is not larger 
than that of the transmitting node after the transmission. This 
pushes the data flow from the source to the destination with 
a positive back pressure. We then define 



7; 



if rec(Z) 7^ e 
otherwise 



Recall that tran(Z) and rec(Z) denote the transmitting 
and receiving node of link I, respectively. Let ei(t) = 
argmax e { 

(*) ~ 9rlc(0,e(*) _ Tf } be me flow on link 

I that has the maximal modified differential backlog, and 



Wl (t) = max [qt n(l)Mt) (t) - <&(,),»,(*)(*) -7r W >0 
nonnegative differential backlog of I at time t. 

The network randomly generates a maximal schedule S(t) 
in a fully distributed manner as in [19], [20], and solve 



is the 



max [ w l(t)(J-l{Pl{t)) - («L(J)(*) + 9r b ec(0 



p(t) ies(t) 
s.t. 0 < 



(tj)W) 



]T Pi(t) < min q b n (t) 
zef2„ue„ 



(ID 



where /x;(P;(f)) = log ^1 + ^p-^j- The objective function is 
separable and concave, and the constraints are local and linear, 
so this optimization problem can be solved in a distributed 
manner [21]. 

Multihop Routing: 
When wi (t) > 0, transmit for flow that is destined to e; (t) with 
rate W (£(*)), i.e., tf {t \P{t)) = W (£(*)) and tf{?{t)) = 
0, Ve^e ; (i). 

It is apparent that MRC and routing can be applied by each 
node independently. 



Observation 1: From the objective function in Equation (11), 
wi{t)m{Pi{t)) is the utility of link I e S(t), when wi(t) 
is large, i.e., the differential backlog of the data queues is 
large, MPA tends to allocate a larger value of Pi(t) so that 
^ttanm ei(t)W w ^ decrease. Further, the routing component is 
used to decrease the length of all data queues over all flows. 
Observation 2: ('? t '' ran ( i )(i) + 9 r f ' ec ( i )(i))-P/(i) is the cost of link I, 
once 9L(I)(*) or is discharged, q* mn{l) (t) or <? r b ec(;) (t) 

tends to increase by Equation (10), then MPA allocates smaller 
Pi(t) so that the discharge event happens infrequently. 
Observation 3: From Figure. 2, as q^. in ^(t) or <Z r b ec( -;)(i) 
increases, (j£an(;)W or ^rec(J)(*) decreases, then MPA tends to 
allocate larger Pi(t) to avoid battery overflow in advance. 
From Observation 2 and Observation 3, MPA expends energy 
neither too conservatively nor too aggressively. 

Note that the amount of actual sensing data R n . e {t) is 
bounded and the endogenous arrivals to any node is also 
bounded due to the peak power constraint and interference. 
Then it is clear that by using the MRC, we make sure 
that the data queues remain within a certain bound and 
this has a positive effect on the batteries as well, since a 
portion of data is not allowed into the nodes. The natural 
question one would ask here is, whether MRC rejects too 
many packets in the first place to synthetically meet the 
constraints. In the following theorem, we show that this 
is not the case. Indeed, if there exists a feasible solution, 
A* = liminfr^oo ^ Y^=o 2~2 n , e #n, e (*) to Problem (A) for 
a given {A(t), t > 0} and {r(t), t > 0}, then the sensing rate 
associated with MRC and MPA can be made asymptotically 
close to A* by appropriately choosing the control parameters V 
and W as the size of batteries B h n grows. We use the notation 
y = O(x) to represent y going to 0 as s goes to 0. 

Theorem 2: If a feasible solution to Problem (A) exists and 
the optimal instantaneous sensing rate vector is R*(t), then the 
joint rate control MRC, power allocation MPA, and multihop 
routing algorithm achieves: 




qle(t)<jV + W, 

q b n (t)<*V + W, 



^ EE *».«(*)} 

t=0 e ) 



(12) 
(13) 



T-l 



> > < lim inf — 



#n,e(*) 



t=0 



K 



-<(2 + /w)-e(W) 



(fv + w-B b n y 

V 



E 

Z6f2„ 



P p 2 eak Lg 
(T 2 



rpeak ' 



w 

O(y). 
(14) 



where a — min/ oi and /i max = log(l 
The detailed proof is in Appendix B. Now, we discuss the 
implications of this theorem. 

In Theorem 2, V is a finite tunable approximation parame- 
ter, which controls the efficiency of the algorithm. W and H 
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are algorithmic parameters, and Bt, is the system parameter. 
Implication 1: W is a parameter in MRC used to bound the 
actual sensing rate. The term e(W) can be arbitrarily small by 
choosing W large enough. 

Implication 2: From Equation (12) and Equation (13), the data 
queues e (t) and virtual queues q h n (t) are all strongly stable. 
Thus, by Theorem 1, p° n < r)°. Given W, as V increases, the 
term O(^) in the gap decreases, but the average data queue 
length and virtual battery queue length tend to increase, which 
implies that the queueing delay may increase. 
Implication 3: Equation (14) compares the performance of our 
algorithm with that of the optimal solution of Problem (A). 
The term J2 n 7 ln(^ + Mmax) captures the influence of network 
disconnection due to battery discharges, and it is small since 
rj° is typically chosen to be very small for all n. 

Implication 4: Given W, the terms 0(— y — and 

O(^) together capture the influence of battery overflow and 
represents the asymptotic behavior of the gap. By choosing V 
and B h n to be large, we can make the performance gap decay 
inversely proportional to the buffer sizes. A larger battery size 
implies smaller amount of wasted replenishing opportunities, 
which is consistent with intuition. 

P LH 

Implication 5: The term — captures the performance loss 
due to imperfect scheduling based on conflict sets. This term 
decreases as H decreases. It is apparent that K < L < oo, 
so the fraction -i is a constant fraction. However, L is a 
very loose bound for K in large networks (large L). From 
Equation (14), we would like both K and H to be small, so 
the question is by decreasing H, whether K will increase to 
a very large value. The following example gives a bound for 
K in a large deterministic uniform network. 

Example: Uniform Extended Network 
In this example, we consider a uniform extended network, 
where nodes are distributed deterministically at a density 
p, homogenously within the network. Let d(m, n) be the 
distance between node m and node n, m,n £ X. Let 
di = d(tran(7),rec(7)), d max = max/gt d u d min = mim e £ d t 
and £ = MH. Let 

(1) 0 < d min < d max < oo, 

(2) there exist a, (3 > 0 such that g m ,n = 
ad _ ^(m,n), Vm,n e N, i.e., signal power decays 
with distance with a path loss exponent j3. 

Theorem 3: For the uniform deterministic extended net- 
work described above, 



K < max { 47rpd max , 6 + £ 



2tt 



arccos 1 



('-*) 



The proof of Theorem 3 can be found in Appendix C. From 
Theorem 3, we can see that K is upper bounded by a constant 
which is not a function of H. The intuition is: as H decreases, 
the conflict area of link I becomes larger and the number of 
links in £/ increases. However, the conflicting opportunities 
of links in £; also increases. Thus, the overall impact of these 
two factors cancel each other and K remains bounded. 



IV. Numerical Examples 

We consider a network topology, shown in Figure 3(a). 
There are 6 nodes, 7 links, and 2 flows with source-destination 
pair (3,1) and (5,2), respectively. We use the rate power 
function as in Equation (1). We choose the power vector of 
the background noise [N] = [2,2,2,2,2,2] x 10~ 14 W, and 
the channel gain matrix 
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and hence a — 0.1. The number of arrivals A% \{t), t > 0 and 
^5.2 (t), t > 0, are modeled as independent Poisson random 
variables with mean A = 20 packets/slot and A max — 30 pack- 
ets/slot. We set ?y°, the threshold of battery outage probability 
to 0.03 for all n £ N and the peak power P pea k = 1 .5 W. The 
backlog threshold 7 = 80 > max„ iE fe rec (0=n W + Ai,e) = 

2/i ^1 + I -^ L ^j, so that the resulting backlog of the receiving 
node is not longer than that of the transmitting node. In all 
simulations, the simulation time is T = 10 6 time slots. 




destination 2 destination 1 

(a) Network Topology 
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(b) Replenishment Process 



Fig. 3. Network topology and a sample replenishing process 

Scenario 1: We first use a replenishment process which is 
formed by a periodic deterministic sine waveform (r max = 1.2 
and period 8000) plus independent Gaussian noise with mean 
0.01, as shown in Figure 3(b) (The cycles imitate the daily 
solar cycles for a solar battery and the average replenishing 
can be simply calculated f = 0.2). All the battery buffer 
sizes are set to be B° — 300J. We simulate three cases: (1) 
crosstalk threshold H = 1 (the notion of H first appear in 
Definition 5), then we obtain K = 3 and there are 7 possible 
schedules; (2) H = 0.05 (H cannot be too small, otherwise 
a link may conflict with all other links. In that case, we may 
need to increase V to a very large value to guarantee a large 
sensing rate, which leads to large queue length and delay 
on the other hand). With this choice, K = 3 (we can see 
that K does not increase as H decreases) and there are 4 
possible schedules; (3) Use the maximal schedules under the 
node exclusive model, then K = 2 and there are 6 possible 
schedules. We choose different values of the control coefficient 
V for the proposed algorithm and compare the results with the 
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optimal value 2 . From Figure 4 (a), we see that as V increases, 
the average total sensing rate keeps increasing and gets closer 
to a value that is much larger than = | (or | for case (3)) 
of the optimal value, which is consistent with Equation (14). 
Case (2) outperforms case (1) since it has smaller H. It is 
not surprising that case (2) also outperforms case (3), since 
case (3) choose schedules according to a protocol model, and 
it results in more schedules than case (2). Time sharing of 
two schedules 2 (these two schedules are contained in all three 
cases) gives the best performance, so the case that has less 
possible schedules should perform better. From Figure 4 (b), 
we see that as V increases, the average data queue length 
(we here only plot the data queue length of node 3 for flow 1 
due to space limitation) keeps increasing but is upper bounded 
by the bound we get in Equation (12). Case (3) has the best 
delay performance, i.e., less average queue length since it has 
smaller K. From Figure 4 (c) we observe that the battery 
discharge probability (we only plot for node 3 here) increases 
to the threshold as V increases. Case (3) goes to the threshold 
slower that (1) and (2) since it has less K which results in 
less virtual battery queue length by Equation (13). Case (2) 
outperforms case (1) in the sensing rate but the discharge ratio 
increases faster, which shows the tradeoff by increasing V. 

Scenario 2: We use a different replenishing processes: r^t) 
and rg(i) are i.i.d binomial random variables Binomial(l, 0.5); 
replenishing at all other nodes are independent binomial 
random variables 0.2 x Binomial(l, 0.5) in even number slots, 
and 0.6 x Binomial(l, 0.5) in odd number slots (F2 = F5 = 0.5 
and f — 0.2 for other nodes), all plus Gaussian noise with 
mean 0.03. We set H = 0.05 and simulate 2 three different 
battery sizes B b = 100J, B b = 10J and B b = 1J (all 
nodes has the same battery sizes). From Figure. 4 (d), we can 
see that the sensing rate increases as battery size increases. 
However, as long as the battery is large compared to the 
average replenishing rate, the improvement diminishes with 
increasing battery sizes. 

V. Conclusion 

In this paper, we studied the problem of energy management 
in rechargeable wireless sensor networks under SINR inter- 
ference model. Our objective was to maximize the average 
data sensing rate subject to QoS constraints on both data and 
battery queues. We developed a mathematical framework that 
allows for multi-hop dynamic routing and handles handle the 
coupling between the data and battery queues. Unlike existing 
works, we allow for a the practical realities of having non- 
stationary replenishment and arrival processes. We develop 
a fully distributed, computationally simple, algorithm that 
achieves a provable constant fraction of the optimum solution 
under the SINR interference model. Through simulations, we 

2 The exact optimal value for Problem (A) is hard to obtain, so we here 
use an upper bound for the optimum that is calculated under the interference 
free environment. For this example, an upper bound for the optimum can be 
obtained by equal time sharing of schedules {1, 4, 7} and {1,5, 6}, and utilize 
the interference free link rate fJ.(f) under infinite buffer and no discharge 
constraint. In Scenario 1, the optimum is 2p(f) = 30.4. In Scenario 2, since 
l 1 (T2) = ^(^5) > X the optimum is A + fJ,(r^) = 35.2 
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Fig. 4. Impact of the control parameter V and threshold H on (a) the average 
total sensing rate, (b) average data queue length of q£ j, and (c) the battery 
discharge probability of node 3 for Scenario 1 . Impact of battery size on (d) 
the average total sensing rate for Scenario 2. 



show the efficacy of our solution. The impact of fairness and 
channel variations (fading model) can be easily accommodated 
into the mathematical formulation. However, an open question 
for future work will be to develop distributed algorithms that 
are provably efficient under a fading model. 
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Appendix 

A. Proof of Theorem 1 

Using the idea similar to [12], we have the fact that if 
any queue represented with Q(t) is strongly stable, then 
limsupy^oo ^P- < 0. Hence, if q\{t) is strongly stable, then 

limsupy^^ < 0. From Equation (10), we have q\{t + 

1) > q b n (t) -v° n + Eien n u e„ ^(*) - r n (t) + M n (t) + l°(t). 
Note that q b n (t + 1) = q b n (t) - £* e n nU e n m) + r„(t) - 
M n (t). By summing from 0 to T — 1, dividing by T and 
taking lim sup of both sides, we have 



lim sup 

T^oo 



> lim 



lim 



ft(0) 
T 

Q b n(0) 



T-l 

Ttf+liminf-X;™ 

T ^°° 1 t=o 

-Q b n(T) 



at slot t + 1. Consider the data queue q d e (t + 1) maintained 
at any node n for flow destined to any node e ^ n at slot 
If node n received data destined to e from other node 
in at slot t, then by the routing policy in Section III and 
definition of w (m ,„)(t), <?l, e (t) - _„,_(*) > 7 ( e m , n) , where 
(m, n) is the link from node m to node n. Choose 7 such 
that the resulting backlog of the receiving node is not longer 
than that of the transmitting node (let /ijjj ax to be the maximum 
endogenous arrivals, then 7 = Ml!!ax+^ satisfy this condition), 
we then have q d e (t + 1) < q^ e (t) + 7, then q d e (t + 1) < 
<7m e (t) < <ZmaxW ^ + W. If node n did not receive any 
data destined to e from other nodes, then it can only have 
exogenous arrivals. Clearly q d n e (t + 1) < q^ e (i) < ^ + W 
if there were no exogenous arrivals. If there were exogenous 
arrivals, by MRC of Section III, we must have q d e (t) < 
then q d e (t + 1) < ^ + W. Thus, Equation (12) holds. 
Proof of Equation (13): Proving Equation (13) from MPA is 
the same as the node exclusive interference model. We here 
prove for a more general case with the rate power function 
fn(P(tj). Consider any link I* e L, when the power assign- 
ment on other links are fixed, fn* (•) = log(-) is monotonically 
increasing, reversible, and differentiable on Pi(t). Hence, we 
have M.(P(t)) < M*(P(t) - Pi-(t)f) + P t *(t) with fixed 
power on other links, where I* is a vector of length L with 1 at 
the T-th index and 0 at others. Since P(t) -Pi- (t)T* € LT(i) is 
the resulting power vector by deactivating link 1* from P(t), 
if P(t) E n(t), P(t) - Pi*{t)T* e n(t) also holds since the 
battery constraints will not be violated. Consider any node 
n* e K and any link I* e Q n , . Since u t , (P(t)) < (P(t) - 
Pi.{t)T*)+Pi.{t) and W (f*(t)) < fii(P(t)- Pi- (t)f*), l?l* 
(for link I ^ I*, energy on I* is its interference, then reduce 
the interference will increase the data rate), we then have 

]T [ Wl {t)vn(P{t)) - (__„,(,)(*) + Pl W 

I 

J2[MtMP(t) - Pi*(t)t*) - 

I 

■Pl(t)] +Wl,{t)Pl*{t). 



< 



?tran(0 



(*) 



rl b 



(*)) 



Further, since P(t) - Pi-(t)l* E U(t), we have from PA 
that 

£ [ Wl (t)m(p(t)) (qi n(i) (t) + _*_(!)(*)) mi 

I 

> ]T [wi(t)»i{P(t) - Pi*(t)T*) - (_*_.(,)(*) + «L.(o(*)) 



0 



0, limsupy^^ SiM 

tr, (0) n 

n j, — 0, so we get 



Since limsupy^^ Hn ^ ' < 
limsupy^^ = 0, and lim T 

p° n = lim sup^^ i EL^o 1 l°n (*) < V° n , Vn e K. ■ 

B. Proof of Theorem 2 

Proof of Equation (12): We prove Equation (12) by induction. 
Let QmaxW be the maximum data queue length for all flows 
at slot t. Assume that q^ ax (t) < + W (holds for t = 0 
by letting <7^. e (0) = 0, Vn, e e 3M)> nee d t0 show that it holds 



Thus, (__„,(,.)(*) + ql^it)) PAt) < Wi,{t)Pi,{t). If all 
the incident links of a node is deactivated, then the virtual 
battery queue this node does not increase anyway. Without loss 
of generality, we only need to consider the case that P;» (t) > 
0. Therefore, $_„_(,.)(*) + &,(,.)(*) < w t .(t) < q d m .^t) < 
can be arbitrary, we have q\{t) < 



KV 

iv 
2 



+ W. Since n* ,1* 
+ W, Vn e N. 
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Proof of Equation (14): Define L(q^(t), ^(t)) = 

~k En,« «e(t)) 2 +wJ2 n (#(*)) 2 ■ From Equation (4) and 
Equation (10), we have 



A(t) = L(^(i + 1), ^(t + 1)) - L(f{t),?{t)) 

^ JC E [<«(*) + 2 Mm a xA„, e (i) + + _! ^ [(1 + 

n,e n 

2 



Proof: Note that log(l + a;) < x for x > 0. Then 



-log 1 + 




01 + P1 + Eiea^i Ki p i ) 



< 



p^f + KY 



+ 2<r„ 



K 



< 



A' 



E [(^(O.e^t) 



<r; + P 

Ppeak 
o 



E 



< 



^E 



hijPi 



hijPi ai + E 



hi.iPi 



E hp^ 



peak Qj 



ieS,ijtl 



(«L.(i)(*) + (*)) + | E [ 7 »W- 

n 

r n (t) + M„(i) . 



Note that if <4 n(i)iei(t) (t) - ^(O.ejwW < *en 
u>z(i) = 0 and /xj(P(t)) = 0 where P(i) is the solution of 
Equation (11), then («L (I) , ei(t) (t)-(4(0,e 1 (t)(*))w(^W) = 
wi{t)m(P{t))\ otherwise, q£ m{l)Mt) {t) - g^. (l)t%(t) (t) > 
wi(t), then (<7 t l n(/)jei(t) (i) - Q Tec (i) !ei ( t ) 
wi(t)[ii(P(t)). Combine with Equation (4) and Equation (10), 
we then have 



the last inequality is from the fact that the link in S has 
crosstalk less than H by the algorithm. ■ 

By Lemma 1, we have wi(t)fj,i(P(t)) - (<7 t fc ran(0 (t) + 



(t))Pl(t) 



- Eles(t) (**(*)) (q b um{l) (t) + 

> Eies(t) [wi{t)n,(W)) - (Cn(Z)W + 

_ EieS(t) (t + Anax)c;. 



Observe Equation (15), it is apparent that MRC is try- 
ing to minimize the term J2 n ,e [i^n.eW - V - ] R n . e {t), 
and MPA is trying to maximize the value of the term 



E 



ies(t) 



By 
e that 



«;,(t)w(fi(t)) - (Cn( 0 W + &c(J)(*))^(*) 

MRC, we have for any n, e G N, n 7^ 

[KeW - ^] Rn,e(t) < [KeW " ^] K,e(*) when 

A n , e (t) < Then, 



A W ^ ^ E K,eW + 2/wA,,e(*) + 2 M 



max 



K 



E 



[(1 + Ppeak) 2 + {V°n) 2 + ^ax + 2<rw] + ^E 



n,e L J n 

~i° n (t) + r„(t) + M„(t)] - I X! [^Ww(^W)- 
(flS».(o(*) + «L(o(*))^(*)] • < 15 > 



< 



< 



R n ,e{t) 

KJt)+ 



K 



K 



<e(*) V 



[A n , e (t)>W] 
^n,e(*) + VA n,e{t)l[A n , e (t)>W], 



Lemma 1: Let S denote the set of activated links chosen 
by the MPA algorithm, then for I € S, we have 



m(P) =log 1 + 



the last inequality can be obtained simply by breaking 
into cases. Since the optimal solution for Problem 
(A) may not be unique, we let V* be the optimal 
solution set and P* € 7* be any optimal solution, for 
Problem (A). In time slot t, let P m (t) be the value 
that maximize the unconstrained objective function 



E 



ies(t) 



Only when P m (t),P*(t) £ 11(f), we may have 



> m{P{) — ci= log 1 



°l + E Je £,^Z ^UP 

PC 



E/es(t) 
E/es(t) 



Q 



«;,(t) W (fi(*)) - (^(t) + ^ ecW (*)) p '(*) 

n(t), 



< 



where c; = Eiesi^; ^^peak is some constant related to 
the threshold h of the crosstalk. 



When P m (i),P*(t) £ U(t), MPA will drain 
at least one battery in order to maximize 

E ;e s (t ) [^)m(to) - (^ m(0 (*)+ tc (0 w) p r(*)" 

within n(£) which is a concave objective function. 
Under this situation, we must have = 1 for some 

n G Ji since r n (t — 1) > 0, Vn e N. Further, note that 
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E 



w(P*(t)) > ni(P*{t)), VI e £, we then have 
A W ^ ^ E K-eW + 2/wA,,e(*) + 2 M Lx] + ^ 

n,e n 
[(1 + Ppeak) 2 + (V°n) 2 + ^ax + + (KV+ 

2A n ,e(*)) (2 + /W)i^(t)] + f£ P„, e (i) + ^ E 

n,e n.e 

|< e (t) v]R* n Jt) - |E E (^an(0 >ei (*)(*) 

d(i), e (t)-7)ft s (^(i))ip 6S(f )] + ^^(t) 



^ P,*(t)-r n (t) + M n (t) + J2 ( V + 

^-4max)ci + V ^2 ^n,e(*) Vn,e(t)>W] 



sides of the inequality 1) > q b *{t) - r„(i) - 

e + 2~2ien ye an( l rearranging terms, we obtain 

[Eien n ue n p i(t)-rn(t)] < q b n (t) [q b n *(t + 1) - 

<7n*(*) + e ■ By summing from 0 to T - 1, dividing by T 
and taking liminfr^oo, we have 



T-l 



E p nt)-r n (t) 

en„\j e„ 



-liminf- V V^(t) 

t=0 n 
T 

: ^ ^ \ E E ^(*)(?t(* - 1) - #(*)) 

^ t=l n 

+ - V liminf gn(r)^*(r)- ^(0)^(0) 



T^oo 
T 



+ liminf-yy 



''max 



l 



- ^ E K,eW + 2(A„, e (i) + 7 ) Mmax + 2 M : 

n, e 

E [(! + P P^) 2 + (^) 2 + ^ax + 2<r max + (^+ 
n 

2A„, e (i))(2 + Mmax )7°(i)] + A^^(t)M n (t)+ 



t=l n 



T 

< ^(1 + Ppeak) liminf + + ^ 



t=l n 



K 



+ E - E )/*f (^(*))^[i e sw] 



+ E (^+ 

ies(t) 



1 , .iVA" AW, 



since the average queue length of the auxiliary queue remains 
finite no matter how large the system parameter B h n becomes 
and is not related to the algorithmic parameters V, H and W. 
By letting e — > 0, we finish the proof. ■ 

Lemma 3: 



1 A )c I 2E " gt(t) 



E *?(*)-»•»(*) 

-ien„u e„ 

+ V5^A n ,e(t)i [Ani .( t)>W] . (16) 



T-l 



V T 



liminf - J2 E ^ W M « W < E °( 



t=0 n 



V 



Lemma 2: 



EE to 



t=0 n 



53 f?(i)-r n (t) 



<0(-). 
Proof: Note that 

j T-l j T-l 

limsup 53 < liminf = k(0) + 53 r„(t) 

t=o ien n u e„ t=o 

T-l 



= liminf^ 51 r «(*)- 

1 — >oo _£ — 



t=0 



Construct an auxiliary queue with the following evolution: 

+ 1) = (*) - r„(t) - £ ) + + 53 i?(t), 

where e > 0 can be arbitrarily small. By Lemma 6, 
q%*(t) is strongly stable. By multiplying q\(t) for both 



Proof: Consider and n £ K. Without loss of generality, let 
9n(()) = 0. Observe Equation (10), we have the following 
cases: 

i) if E/ ei2 „ue„ P 'W > r n (t), l°(t) - 0 and q b n (t) > 0, 
then M n (t) = 0, q b n (t + 1) < q b n (t) and (t + 1) - g*(t) < 
qn(t) — q^(t + 1), i.e., even if q\(t) increases, the increment 
is no larger than the decrement of <j£(t); 

ii) if Ei £f2 „ue„ P M < r n(t), l° n {t) - 0, then if q b n (t + 
1) > 0, q b n (t) - q b n {t + 1) = r n (t) - E ie o nU e n ^(*) " 
M n (t) - V ° n + « - gt(t))+ > ^(t + 1) - ( t ) = r n (t) - 
E/Gf2„U e„ p l(t) ~~ M n (t), i.e., the decrement of q\(t) is no 
less than the increment of q b n {t), else if q\(t + 1) = 0, it goes 
to case iv); 

iii) if l° n (t) = 1, then q£(t+l) = min[r„(t), S*] by definition 
of discharging event Equation (3). If M n (t) > 0, then this 
means r„(t) > P^ and -r„(t) + M n (t) = -B b ; 

iv) if = 0, then q b n (t)M n (t) < 0 anyway. 

Combine the above discussion with Equation (13), we have 
that Vt > 0, if M n (t) > 0 and q b n {t) > 0, we must have 
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q b n (t) < + W) + max[r max , 1] - B b n . Thus, 



T-l 



t=0 



< 



(( 



KV 



W) + maxlrw, 1] - B b n )+r n 



V 



0{ 



(fV+W-B b n )+ 



V 



)• 



Under the optimal policy, the amount of data destined to the 
same destination may be routed through different routes. We 
first logically decompose the flows that are differentiated by 
destination into flows that are differentiated by routes under 
this optimal policy. Let 2f denote the set of flows after the flow 
decomposition, then any flow / e J has a fixed route. Let e(/) 
denote the original flow from which flow / is decomposed, 
and Uf(l) denote the one hop upstream link of link I in the 
route of flow /. Let H* = {H l j*} denote the routing matrix 
under the optimal policy, where Hj* = 1 means the data of 
flow / is routed through link I. We further logically decompose 
the data queues at nodes into data queues at links under this 
optimal policy. Since the sum queue length of the decomposed 
queues equals the original queue length, we have 



\iee„ ien„ 



— £ ^ttan(;),e(/)(*)^ttan(/),/(*) ^ (q£m(1) ,e(/) (*) ~ 
U L 

d(0,e(/)(*))A*f(^(*))^H6S(t)]' 

_V^rf m \Y,n R *nA t ) Hl f* 

~ 2^ 9 tran(;),e(/)W 

U 



-n{(P*(t))i [lesm - 



r* (t)H^ f ^* 

; £ ^(o,e ( /)(*) (Qfj(t + 1) - oft®) 

u 



W(0^ Wi/ (i + l)-Q^( 0i/ W) (17) 



where (5 > 0 and Qf^t), VZ £ £, / € J has the following 
evolution: 



Qft(t + 1) = (<3fr(f) - /if (P*(t))7 [l6S(t )] -S) + 

+ ^£k,/(*)4*- 



(18) 



Lemma 4: If all the queues Qfj(£), VZ e £,/ e 
3" with evolution Equation (18) are strongly stable, 



then ^ liminf ^ £ t=0 E„, e ?n,e(*) 
-£* e n>f(^*W)^es«]l <0(f) 



iee„ 



Proof: Since < < e (i) - e (t + 1) < Mn 
0, Vn, eeN, then 



Vt > 



T-l 



liminf ^ 

T^oo T 



U 



t=o 



~ ^i(l),e(/)(*) (Qu f (l)J*( t + 1) _ QlJfWJ*^)) 

u 



< 



E 



lim 

T^oo 



■*tran(/),e(/) 



(T)Qt,f*(T) 



+ /J, m - dx liminf 



T^oo T 



1 T 



t=i 



1 T 



+ lim 

T^oo 



^■(i), e (/)(0)^ /( i),/.(0) 



T 

T 



< E k« h T ^ f f £ of,/. (*) + w ^ £ 



*=i 



T^oo T 



t=l 



Qu f (i)jM 



(19) 



The average queue length of Qfj(t) is not related to V, H 
and VK. Combine with Equation (17), we obtain the result. ■ 
Now we only need to show that all the accumulative data 
queues Qfj, VZ G £*,/ G J are strongly stable. We know 
that all the data queues are strongly stable under the optimal 
policy queue evolution, then the queues are still strongly stable 
after flow and link decomposition. This further implies that 
the accumulative queues are strongly stable (note that the 
accumulative queues being stable does not imply each queue 
in the routes are stable), i.e., the queues with the following 
evolutions are strongly stable: 

Qt f (t + i) = (Qf; f (t) - »{(P*(t))) + + e</(*)4*- 

n 

(20) 

Let Ji,{* = liminfT^oo ^ Y,J=o A and X i* = 

liminfr^oo \ J2t=o E n ^n,/(*)-^/*- Without l° ss °f gen- 
erality, assume p,f* > 0. We know the optimal sensing rate is 
feasible, so Xf* < fif*. Then for the normalized queues with 

the following evolutions, the arrival rate is also feasible: 

m 



V /* 



+i) = (<&(*) + i 



(21) 



Note that l[i£S(t)] i s a random variable which is independent 
of P*(t). Thus, l[i e s(t)] and (P*(t)) are independent and 
then uncorrected, so 



T-l 



T-l 



; ,,..), = jif* liminf I^7 [;e ,, : , . 

t=o 

(22) 



t=o 



In order to show that all the queues with evolution Equa- 
tion (18) are strongly stable, we need to show that the 
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normalized data queues under the MPA algorithm with the 
following queue evolutions are strongly stable: 



1 



[ies(t))j 



f* 



f* ■ 



(23) 



Lemma 5: If the arrival rate is feasible for queues 
with evolution Equation (21), then the queues with evolution 
Equation (23) are all strongly stable. 

Proof: Similar to the idea of [20], we let L(q d *(t)) = 

E U Qtf( t )(Ekei[J^J& t ))- Let £ < denote the set of 
links that will not be activated under MPA algorithm when 

link I is activated. Then 



A L ^, ) (t)=L(f*(t + l))-L(q d *(t)) 

<2 e m\E e / A,,i 



u--qf;,(t)>i 



1 feG/USi 



-i 



[fces(t)] 



L 2 F 2 I 1 



max 




< K 



Note that EkeiU e, hkesd)] > 1 and EkeiU \ e, ^ 
for any / G 9\ By summing from 0 to T — 1, dividing by T 
and taking liminf 7^ co for the above equation, we have 



limsupi]T£X/(*) 

T ^°° t=0 ij 

1 T_1 

<l + limsup-E E 



T^oo 



00. 



t=0 U:qf;,(t)>l 



By Lemma 5, queues with evolution Equation (23) are all 



strongly stable, then < liminfr^oo \ J^Uo 7 [ies(t)]- 

Combine with Equation (22), we have 

^^p^eY^e^/^/I 



T-l 



< h T ^ f ^ E (^(^*(*)) j [ies(t)] + • (24) 
t=o 

The following lemma gives the sufficient condition for 
strong stability under sample mean argument. 

Lemma 6: For any queue with evolution q(t + 1) = 

(q{t) - + A(t), if A(t) < A max , Vi > 0 

and limsupy^^ ^ Xt^) 1 X ( l ) < liminf t^oo y Et^) 1 /"(*)> 
then is strongly stable, i.e., limsup T _ >00 ^ X^^Jo 1 < 
oo. 

Proof: For any time t, either there exists a time J(i) < t 
such that g(i + 1) = q(i) - n(i) + X(i), Vi € or 
«(t+l) = A(t). Further, limsup^^ ± Ej=o (W) ~ < 
0 implies that there exists T 0 such that Y^t=T x ( A W — M*)) < 



oo, V7i,T 2 > T 0 . Then, 

1 T_1 
limsup- E «(*) 



T^oo 



t=0 



1 T_1 

< limsup - V I A 



t 

max "F ^ ' 

i=J(t) 
T-l t 



(A(i)-M^ 



< A max + lim sup ^ ^ 



(A(i)-MO 



< oo. 



Let a = min; ct/, then Xies(t) c £ — ^ Pp ^ kL ^ t7. Combine 
with Lemma 2, Lemma 3, Lemma 4, Lemma 6, by summing 
from 0 to T — 1, dividing by T and V, taking liminf for 
Equation (16), we get 

n L t=0 e J 

> V J luninf ^ E E ^ ~ < C + *») " 

t=0 e 



'( 



O' v2 



V 



^peak-^ 



T-l 



-E h > inf ^ E A n,eO)^„, e (t)>W/]. 
— i — >oo 1 — 

n,e t— 0 

Since liminf T ^oo ^ Et^) 1 Ai,e(*) < ^max < oo and 

^ T-l ^ T-l 

^inf t E - liminf - E A n , e {t)l [A {t) < w] 

t=0 
1 

+ liminf - E A n,e{t)l [An<e {t)>w]- 



t=0 



T^oo T 

t=0 

By Monotone Convergence Theorem, we have 
liminfr^oo \ J2j=o ^«,e(*) 7 [A„, c (t)>w] —> 0 as W -> oo, 
which implies Equation (14). ■ 

C. Proof o/ Theorem 3 

Consider any link I e £, for Vfc G we have /ifc^ = 

ginm(fc),, C c(Q _ ( rf(tran(fc),rec(Q) \ 0 > jj 0 T hi k = 9 "'"'<''' rec <'°> = 
^(tran(0,rec(fc)) ^-' 3 > ^ i g d(tran(fc ) ; rec(/) ) < dl H~^, OX 

d(tran(Z),rec(fe)) < d k H~K 

As shown in Figure. 5, any link whose transmitter is in the 
diH~e neighborhood of rec(Z), and any link k whose receiver 
is within the distance dkH~P of tran(7) are conflicting with 
link I. We need to find an upper bound for the number of 
nonconflicting links in Let r; = diH~?, V7 e L, r = 
d m inH" ? and R = d m . ix H~ ? . i 

I) We first consider the diH~P neighborhood of rec(7). 
Fact 1: From Figure. 5, it is apparent that with a fixed H, if 
tran(i) and tran(j) are both within as d(tran(i),rec(/)) and 
c?(tran(j), rec(Z)) increase, the angle 9 gets smaller by placing 
tran(z) and tran(j) closer without letting i conflict with j. 



Thus, we can just assume d[ = d max and 
Fact 2: If rv < cL and r, < cL, i.e., H > 1, the transmitters in 
the c?;iJ 0 neighborhood of rec(Z) can be placed very dense 
and the upper bound for the number of nonconflicting links in 
this area is 2tt pd"^. dx H~ ^ < 2k pd 2 ^. 

We then consider the case when H < 1, from Figure. 5, it 
is easy to have the following 



(d, + R- Rcos 0 R ) 2 + R 2 sin 2 9 R = r 



2 



Rdj 

^max 



Then cosOn = 1 - R ^ + _ i), and it is 

minimized by letting di = d m i n . Thus, 



COS 0ft = 1 



j2 r>2 _ 12 

2d 2 m , x R(d min + R) 



d 2 - d 

"mm I ^ 



2rf max I dmin-ff 3 ^min + dmax# 

which is monotonically decreasing to 1 — j^f- as H — > 0. 

Thus, 0 m j n = arccos ( 1 — ^Jf- ) and there are at most 
transmitters can be placed along the outmost circle of rec(7)'s 
neighborhood. By Fact 1, in the inner circles, the radius 
gets smaller, and the minimum angle for placement increases. 
Therefore, I 1 is also an upper bound for number of 

"min 

transmitters placed along all the inner circles, then there are at 



most £ 



nonconflicting links in the diH 



neighborhood of rec(Z) when H < 1. 

II) We now consider the conflicting links induced by tran(7). 
When H > 1, Fact 2 is also valid here. So we consider H < 1. 
From Figure. 5, it is easy to see that placing a receiver of a 
link with length d m . ix and a receiver of a link with length 
d m i n alternatively makes the angle larger. So, we only need 
to consider placing all receivers of link with length d m . dx or 
placing all receivers of link with length d min . Further, the 
conflicting distance to tran(Z) is determined by the length of 
other links, so we only need to take care of one circle. Similar 

to the calculation in I), we have cos ip R = \ + ^ = 1+ 2^ = 
i _|_ 4a _ cos p r which is monotonically decreasing to \ as 
H — > 0. Thus, 9 m i n = | and there are at most 6 receivers can 
be placed around tran(Z). 

Combining I) and II), we finish the proof. ■ 



